Numerical Calculation of the Functional renormalization group fixed-point 

functions at the depinning transition 
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We compute numerically the sequence of successive pinned configurations of an elastic line pulled 
quasi-statically by a spring in a random bond (RB) and random field (RF) potential. Measuring 
the fluctuations of the center of mass of the line allows to obtain the functional renormalization 
group (FRG) functions at the depinning transition. The universal form of the second cumulant 
A(u) is found to have a linear cusp at the origin, to be identical for RB and RF, different from the 
statics, and in good agreement with 2-loop FRG. The cusp is due to avalanches, which we visualize. 
Avalanches also produce a cusp in the third cumulant, whose universal form is obtained, as predicted 
by FRG. 



Universality is often more difficult to characterize in 
random systems than in their pure counterparts. Sample- 
to-sample fluctuations complicate the analysis and the 
nature of the critical theory may be different. One promi- 
nent example is the zero-temperature (T = 0) depinning 
transition from a pinned to a moving state, which occurs 
when an interface is pulled through a random medium 
by an external force / beyond a threshold f c . Its un- 
derstanding is important for magnets ferroelectrics 
3|, super-conductors 0, 0, density waves 0, wetting 
aj, dislocation Hand crack propagation ||, and earth- 
quake dynamics [2J- At the transition the interface dis- 
placement u(x) is expected to scale as u(x) — u(0) ~ x^, 
where x is the d-dimensional internal coordinate and £ 
the roughness exponent. The analogy with critical phe- 
nomena, suggested by mean-field theory [lol |. was ana- 
lyzed using the functional renormalization group (FRG) 
to one loop [flj, [T2J ■ 2-loop FRG studies resolved the ap- 
parent contradiction that statics (/ = 0) and depin ning 
(/ = fc) cannot be dis ting uished at one loop naiij. 
In the earlier works [1(1 Il2j , the presence of a diverging 
length scale was argued to lead to the universal behav- 
ior observed at the transition. This correlation length 
was observed in numerics for the steady-state d yna mics 
above f c , but only in transients below f c [lEl Ho . 1171 Il8| . 
The FRG study |lj| of thermally activated motion below 
threshold showed a more complex picture with additional 
length scales involving both statics and depinning. This 
is in agreement with a recent numerical analysis of the 
T = + steady state in that regime; |2fJ shows that (i) 
there are no geometric diverging length scales at /~ for 
this steady state and (ii) the roughness is given by the 
equilibrium static exponent at small scales and by the de- 
pinning exponent at large scales for all < / < f c . The 
physics is thus more subtle than in standard critical phe- 
nomena. The 2-loop FRG is a good candidate to describe 
this physics as it contains a mechanism for a crossover be- 



tween statics and depinning directly in the quasi-static 
limit T, v — > (by the generation of an anomalous term 
in the (3- function at any / > 0). It is thus important 
to directly test the central ingredients and properties of 
the FRG approach in the dynamics, making contact with 
observables beyond critical exponents. 

Recently a method to measure the fixed-point func- 
tion of the FRG for the statics of pinned manifolds 
was proposed |2l|. Exact numerical determination of 
ground states for interfaces in various types of disor- 
ders [2^] shows a remarkable agreement in the statics 
between the measured renormalized pinning-force corre- 
lator, A (it), and the 1- and 2-loop predictions from the 
FRG [ii III Q. This method has been extended to 
the quasi-static depinning [2j|. The aim of the present 
paper is to compute numerically these fixed-point func- 
tions for depinning. Outstanding predictions of the FRG 
which we test here are: The existence of a linear cusp for 
A(tt), a single universality class for both random-bond 
(RB) and random-field (RF) disorder, the difference of 
A(tt) between the static and depinning fixed points, and 
a comparison with 1- and 2-loop predictions. In addition 
we study the third cumulant, which also exhibits a cusp. 
The cusps in these FRG fixed-point functions can directly 
be related to "avalanches" or "dynamical shocks" . 

The main idea to study depinning, described in |2fij . is 
to put the system in a quadratic potential and to move 
its center, denoted w, monotonously and quasi-statically: 
The difference between the center of mass of the mani- 
fold and w will fluctuate, and its second cumulant yields 
precisely the function A (it) defined and computed in the 
FRG. In the continuum, the zero-temperature Langevin 
dynamics is described by the equation of motion: 



d t u(x,t) = T w u){x,u{x,t)) 



(1) 



T w {x, u(x)) = m 2 (w — u(x)) + cV 2 u(x) + F(x, u(x)), 
where !F w {x, u) is the total force acting on the man- 
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FIG. 1: Structure factor of the line (L = 512) for RF disorder 
(curves are shifted for clarity). The crossover between the 
depinning slope and the plateau determines the correlation 
length (oc m _1 ). For m — 0.01 the correlation length is bigger 
than the system size. 




ifold, c is the elastic constant and 



is the curva- 



ture of the quadratic potential which acts as a mass for 
the field u. F(x,u) is the random pinning force. For 
RF disorder, F(x,u) is short-ranged with correlations 
F(x,0)F(x',u) = A Q (u)S d (x -x'). For RB disorder this 
force is derived from a short-ranged random potential 
V(x,u), F{x,u) = —d u V{x,u). 

Starting from an arbitrary initial condition u in ; t (x), 
and a given w = u>o, the manifold moves to a locally sta- 
ble state u Wo (x), i.e. a zero-force state J r Wo (x, u Wo (x)) = 
0, which is stable to small deformations. Increasing 
w, u(x) increases slightly (and smoothly if F(u,x) is 
smooth in u), while the configuration remains stable. 
At some w = wi, the state becomes unstable and the 
manifold moves until it is blocked again in a new lo- 
cally stable state u Wl (x). We are interested in the cen- 
ter of mass (i.e. translationnally averaged) displacement 
u(w) = L~ d J d d x u w (x). The function u(w) exhibits 
jumps at a discrete set of values of w and is in gen- 
eral dependent on the initial condition. However, due 
to the no-passing rule |26j, we can prove that there ex- 
ists a w* > wo such that the orbit u w>w *(x) becomes 
independent of the initial condition Ui n n(x), and wq. A 
stationary state is thus reached after a finite w — Wo, on 
which we focus. 

We check these predictions numerically for a string 
(d = 1). To solve Eq. QJ, we discretize the string along 
the x direction, x — > j = 0, . . . , L — 1, keeping u „,{j) 
as a continuous variable. A very efficient algorithm |2j| 
finds the exact location of the succession of locally sta- 
ble states. For RB disorder, we generate, for each j, 
a cubic spline V(j,u(j)) interpolating a large number 
(10 2 — 10 3 ) of uncorrelated normal random points, of reg- 
ular spacing a, zero derivative being imposed at the first 
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FIG. 2: Finite-size study of f c (m). The extrapolation f c = 
/ c (0) corresponds to the critical force for the infinite system. 



and the last points. Once u(j) passes the last point, a 
new spline is generated. For RF disorder, F(j,u(j)) is 
taken as a linear interpolation of regularly spaced nor- 
mal random points. The discretization in x, in the limit 
a — > 0, preserves the statistical tilt symmetry (STS) of 
the continuum model. (Only very small corrections to 
c are expected as ma/^/c <C 1). The Fourier modes 
and center of mass of the discretized line are defined as 
u q = Ylf^o e lq] u w (j) and u(w) = u /L. 

We have observed that in the transient regime w—u(w) 
increases on average linearly with w and reaches a plateau 
in the stationary state. There, the line is depinning-like 
rough up to a scale of order 1/m where the confinement 
due to the mass takes over. This is apparent on the 



disorder-averaged structure factor S(q) 



plotted 



in Fig.njfor various masses: it exhibits a plateau at small 
q. In the steady state the fluctuations of w — u(w) are 
related to the FRG functions. The first cumulants are 
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FIG. 3: Universal scaling form Y(z) for A(u) for RB and RF 
disorder. 
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FIG. 4: The difference between the normalized correlator 
Y(z) and the 1-loop prediction Y\{z). Averages over 10 7 - 
10 8 samples were performed. 



FIG. 5: Data collapse for the universal function Q(x) defined 
in 10, for RB and RF disorder. The line represents the 1-loop 
prediction Q(x) = (1 — x) 2 . 




m 2 [w — u(w)] — f c (m) (2) 
to 4 [w — u(w)][w' — u(w')] = L~ d A m (w — w'). 

Since the correlations of w — u(w) decay over a finite scale 
in w, the disorder averages in can be determined as 
translational averages over w. A prediction of the FRG 
is that in the limit mL — > oo the quantities f c (m) and 
A m (w) in (J2J become L-independent. Here Fig. ^ shows 
that this holds for Lm > 5-s/c, as one can check that 
several modes are confined and the correlation length is 
smaller than the system size. The FRG also predicts 
that: 

fc(m) =f c + Cl m 2 -< (3) 
A(u) = m e ~ 2< A(um c ) (4) 

where A (id) goes to a fixed point as m — > (e = 4 — d; 
here d — 1). 

We have studied the behavior of the critical force 
/ c (m) for the two classes of disorder. From 10} one has 
y A(0)m ~ m 2- ^, yielding a parameter- free linear scal- 
ing shown in Fig. [3 For large m the scaling is non- linear, 
while for smaller m it is linear up to the scale where the 
correlation length becomes of the order of L (mL w 5). 
The critical force of the infinite system is defined here in 
an unambiguous way, as f c = f c (m = 0). The resulting 
ci < 0; hence the average force exerted on the manifold 
is smaller than f c . One can see on Fig. El that the slope of 
the two curves coincides. This is consistent with the FRG 
which predicts that it is a universal amplitude, depend- 
ing on microscopic details only through the renormalized 
elastic coefficient c; here c ~ 1 for both models of Fig. 
The study of this and related amplitudes is deferred to a 
future publication. Here we focus on parameter-free fully 
universal functions. 



We now turn to the determination of the fixed-point 
function. Since there are two scales in A(u), we write 

A(u) = A(o)y(u/« c ), (5) 

where Y(0) = 1 and one determines such that 
f Q dzY(z) = 1. The function Y(z) is universal and 
depends only on the dimension of space. We have de- 
termined Y(z) from our numerical data both for RF 
and RB disorder. For small masses the two functions 
coincide within statistical errors. This is visualized in 
Fig. The prediction from the FRG is that Y(z) = 
Y 1 (z) + eY 2 {z) + 0{e 2 ) with e = 4-d. The 1-loop function 
is the same as for the statics and given by the solution 
of jz = y/Yi — 1 — In Yi and 7 = f dy \/y — In y — 1 w 
0.5482228893. Since the measured Y(z) is numerically 
close to Yi(z), as in the statics, we plot in Fig. 0] the 
difference Y — Y\, The overall shape of the difference 
is very similar to the one obtained for the RF statics in 
d = 3, 2, 0, which exhibits only a weak dependence on 
d. However the overall amplitude is larger by a factor 
of about 1.25, both in the numerics and in the 2-loop 
FRG. We have plotted the function Y 2 {z) = £Y(z)\ e=0 
which, as for the statics, is close to the numerical result. 
We also observe a cusp, i.e. Y'(0) = -0.816 ± 0.004 for 
RF and F'(0) = -0.815 ± 0.005 for RB. FRG predicts 
y'(0) = -0.775304- 0.0412061e+ 0(e 2 ). 

To investigate deeper the validity of FRG, we measure 
the third cumulant 29], defined as 

m 6 (w' - u(w') -(w- u(w)))3 C = L- 2d S(w' - w) . (6) 

The lowest order prediction [2l| is S(w) = 
^A'(w)(A(0) - A(w)). To check the scaling in a 
parameter-free way, we define 

i> w poo 

Q(AH/A(0)) := J S(w')dw'/ J S{w')dw'. (7) 
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FIG. 6: Shocks (jumps) in the evolution of the center of mass. 
Decreasing m, several small shocks merge into few bigger ones. 



Q(x) is expected to be universal. Indeed we find, that RB 
and RF give, within statistical errors, identical results, 
close to the 1-loop prediction, see Fig. 

It is instructive to visualize the shape of the function 
u(w) in a single environment as a function of the mass to, 
as shown in Fig. [5] The analogous function in the stat- 
ics, i.e. the position of the center of mass of the ground 
state, exhibits shocks. In d = the evolution of these 
shocks as to is decreased follows a ballistic aggregation 
process described by the Burgers equation. The "dy- 
namical shocks" or avalanches also follow an interesting 
dynamics which remains to be studied in detail. The fact 
that they do not accumulate, apparent on Fig. |SJ is con- 
sistent with the linear cusp found in the second and third 
cumulants. 

To conclude we have analyzed the dynamics of a man- 
ifold at the depinning transition in a geometry which 
allows a precise and unambiguous comparison to the 
predictions from the Functional RG. By moving the 
quadratic well quasi-statically, we cleanly define the 
avalanches at the threshold. The center-of-mass fluctua- 
tions become universal and are described by the FRG 
fixed-point functions. The main and non-trivial pre- 
dictions of 2-loop FRG are confirmed, namely a scale- 
invariant fixed-point function A(u) with a linear cusp and 
a single universality class for RB and RF disorder. The 
precision of the data allows for a quantitative comparison 
with the statics, in agreement with FRG. A more detailed 
analysis of other universal observables is deferred to a fu- 
ture publication [2^] . We hope that this study opens the 
way to more precise comparisons between theory and ex- 
periments. 
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